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Abstract 

We determine the small signal gain and noise response of an amplifier based on the nonlinear 
response of a quantum nanomechanical resonator. The resonator is biased in the nonlinear regime 
by a strong harmonic bias force and we determine the response to a small additional driving signal 
detuned with respect to the bias force. 
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I. INTRODUCTION 



Almog et al.[T] recently demonstrated classical noise squeezing and amplification in a 
nanomechanical resonator with a significant fourth order nonlinearity in the elastic poten- 
tial energy. This system is dynamically equivalent to the Duffing oscillator, which exhibits 
a fixed point bifurcation: a loss of stability in a fixed point as a control parameter (in this 
case a driving force) is varied[2]. Near the bifurcation, the system becomes very sensitive 
to fluctuations, a sensitivity that can be harnessed to make a signal amplifier [3], or sensi- 
tive charge detector |4j. Such a bifurcation amplifier has recently been used in the study of 
superconducting qubits[S]. In this paper we calculate the small signal gain and noise power 
spectrum for quantised quartic nonlinear nanomechanical oscillator when the harmonic driv- 
ing force is subject to a perturbation. Our treatment complements the study of Buks and 
Yurke^ in which the response of a nonlinear resonator to changes in the linear frequency, 
due to a mass perturbation, enable a sensitive mass detection transducer. 

The Hamiltonian for the nanomechanical system is [7] 



where m* is the effective mass of the nanomechanical resonator, ujq is the linear resonator 
frequency taking into account the applied strain. Kozinsky et al.[2j give the nonlinear 
parameter as 



where Q is the quality factor of the resonator and Oc is the critical amplitude at which the 
resonance amplitude has an infinite slope as a function of the driving frequency. 

The position and momentum operators of the nanomechanical resonator may be written 
in terms of raising and lowering operators, a\ a, 



We also include a harmonic driving field at frequency uip, and amplitude e^, which we call 
the pump, that is used to set the operating conditions of the device, and a weaker harmonic 
signal driving field at frequency Ug and amplitude e^. The Hamiltonian may then be written 
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as 

H = hujQO^a + h—{a + a)Y + h{2ep cos{ujpt) + 2es cos{uJst))(a + a^) (5) 

where x gives the nonhnear dispersion. It is given in terms of the nonhnearity parameter a 
as 

In the example of a doubly clamped platinum beam in [2] , x ~ 3.4 x 10~^ s~^. 

Moving to an interaction picture at the pump frequency, Up, and assuming that ujq, Up » 
X to neglect rapidly oscillating terms in the quartic term, the Hamiltonian becomes 

Hj = hAa^a + hx{a^fa^ + htp{a + a)) + hts{ae'^' + a)e-'^') (7) 

where A = c^o — ^p-, 5 = Us — ojp are the detuning of the resonator and the signal from the 
pump respectively. 

The actual physical model for dissipation in quantum NEMS is still under investiga- 
tion. For the purposes of this paper we will model damping via the quantum optics master 
equation [9j. This assumes that the nanomechanical resonator is under-damped due to a 
weak coupling to a bath of harmonic oscillators. It further assumes the validity of the rotat- 
ing wave approximation for the interaction between the NEMS degree of freedom and the 
bath oscillators which is expected to be a good approximation for sufficiently high NEMS 
frequency and sufficiently weak coupling to the bath. 

The master equation in the interaction picture for the system is then given by 

= — — [Hi,p\ + — (n + l){2apa) — a)ap — pa) a) + — (2a^pa — aa) p — paa)) (8) 

where 7 is the rate of energy loss from the resonator and n is the mean phonon number in 
a bath oscillator at frequency uoq. We will henceforth assume low temperature operation so 
that n ^ 0. In the case of no signal field, this model was introduced long ago in quantum 
optics to describe optical bistability due to a Kerr nonlinear medium |10]. The system has 
a steady state, or fixed point, which can change stability as the driving field is varied. It is 
this dependance of fixed point stability on driving field that can be used to amplify a weak 
driving signal. 

In addition to this dissipative channel, the transducer for the nanomechanical displace- 
ment itself provides an open channel by which the resonator is coupled to the exter- 
nal world. A variety of transducers are possible, including single electron transistors [12], 



super-conducting single electron transistors [13] and super- conducting co-planar microwave 
cavities[T3]. We do not model a specific device here but proceed with a more generic de- 
scription. 

We will assume that the transducer is described by a multi mode bosonic field, with 
positive frequency components, ain{t) , aout{t) , and a carrier frequency close to the pump 
field frequency of the resonator. If we further assume that this coupling is linear and make 
the rotating wave and Markov approximation, the bosonic field after the interaction, aout{t), 
with the resonator is related to the complex amplitude of the resonator and the input field 
to the transducer aj„(t) by [5] 



where aj„(t) is the quantum noise operator for the input noise to the transducer field, which 
for simplicity we will take to be in a coherent state with a monochromatic component 
at carrier frequency uJs and amplitude e^. This assumes that the carrier frequency of the 
transducer field is large compared to ksT/h. We have in mind a transducer in the form 
of a microwave field in a superconducting co-planar strip line [15]. At a carrier frequencies 
of GHz and milliKelvin temperatures such a transducer may be regarded as in a coherent 
state. Under these assumptions we see that the transducer field is simply proportional to 
the operator a{t) for the nanomechanical resonator. The total damping rate, 7 should be 
regarded as including the transducer damping rate 7t. 

In order to see the displacement of the resonator we need to make a phase sensitive 
measurement of the transducer field. This requires a frequency and phase reference. In the 
case of electromagnetic fields this would be done using homodyne detection of a quadrature 
phase amplitude of the field. We are thus led to define the quadrature phase amplitude 
operators for a nanomechanical resonator in the Heisenberg picture as 



where 6 and ujio are the phase and frequency reference provided by some form of local 
oscillator. In our case, the input signal to the amplifier is set to uJs- As we will see, this 
leads to a long time solution in which the signal has Fourier components at frequencies 
ujp ± 5. With this in mind we then assume that the local oscillator frequency is set so 
that ujlo = ± (5. The operators a{t),a'^{t) appearing in these expressions are multi-mode 




(9) 



Xeit) = a(t)e'(^+"'°*) + at(t)e-*(^+"'°*) 



(10) 
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operators describing the quantum noise of the detected signal propagating away from the 
local system. It is conventional to write them in terms of the positive frequency components 
a{u!) as 

/•oo 

a{t) = / duja{uj)e-''^^ (11) 







where [a{u!), a^{uj')] — 5{u! — uj'). The quadrature operators Xe{t) and A'5)+^/2(t) are canon- 
ically conjugate multimode operators with equal-time commutation relations, 

[Xe{t),Xe+^/',{t)\^2i (12) 

We are working in an interaction picture defined with respect to the pump frequency ujp so 
we can write 

Xe{t) = ai{t)e'^^^^'^ + a\{t)e-'^^^^'^ (13) 

where ai{t) — a(i)e*'^p* and we have defined the local oscillator frequency as culq — ujp± 5. 
In terms of the positive frequency components the interaction picture amplitude may be 
written as 

/•oo 

ai{t) = / duja{uj)e-'^''-''''^' (14) 
Jo 

We expect that the states of interest are such that n(a;) — (a^(a;)a(a;)) is significantly 
different from zero only over a bandwidth B centered on a; = a;^, with cUp » B. In an 
experiment B is the order of the resonator line width which is of order MHz, and cUp is at 
GHz. We thus make a change of variable e — cu — cUp, and setting a;^ — > — oo we can write 

/oo 
ded{e)e-''* (15) 
-oo 

indicating that a(e) = a{u!p + e) is the Fourier transform of the interaction picture operator 
ar{t). In a similar way we find that 

/oo 
ded\e)e-''' (16) 
-oo 

but with the definition a^(e) = a{u!p — e). Note that the commutation relations for the 
Fourier components is not the same as that for the original positive and negative frequency 
components, 

[d{e),d\e')] = 6{e + /) (17) 
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II. BISTABILITY 



We first establish the operating conditions for the device when there is no signal present, 
= 0. Under certain conditions, the energy in the nano mechanical resonator as a function 
of the driving intensity can exhibit multiple stable fixed points and hysteresis. We will follow 
the approach of Drummond and Walls [lOj which is based on the positive P representation 
of the density operator. In the semiclassical approximation, the equation of motion for the 
mean amplitude, a = (a), is given by 



a 



-iep- (7/2 + i(A + 2x|an)a (18) 



where the dot indicates a time derivative. The fixed point (or semiclassical steady state) is 
defined by d = 0, which corresponds to a complex amplitude ao must satisfy 



Jp = no 



(19) 



^ , (A + 2xno 



where Ip = ei is proportional to the pump power driving the nanomechanical resonator and 



p 

|2 



^0 = Ic^ol determines the average energy in the nanomechanical resonator hy E = fiw^n. 
Considered as a function of no, Ip is a cubic with turning points at the values of no that 
satisfy 

^ = ^ + ( A + 6xno) (A + 2xno) = (20) 



However when we regard no as a function of the pump power, it is multi- valued and Eq.(20) 
defines values at which the slope diverges, indicative of a change in stability. 

In Figure [T] we plot no versus the pump intensity for various values of A. Clearly under 
some conditions no becomes a multi valued function of ep. In fact it can be shown that this 
will occur for negative detuning, A < 0. Not all the fixed point solutions are stable. To 
determine stability we linearise the equations of motion around the fixed points by writing 
a{t) = tto + 5a{t). The equations of motion for the fluctuation field 5a{t) are then given by 

d'^ \ 5a* j \ 6a* j 

where 

/ -2-i(A + 4xno) -iC , 

M = \ ^ ^ ^ (22) 
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(b)x=0.1, A= V3 




where G = 2x^0 and ao is the solution to 

"7 



Then we can write ao = \/noe^'^° where 



2 + 2(A + 2xno) = -tep (23) 



tan (f)o = (24) 

As we have taken ep as real this is the phase shift of the oscillator field from the pump field. 
The eigenvalues of the linearised motion determine stability. These are given by 

= -| ± V(A + 6xT2o)(A + 2xno) (25) 

For stability the real parts of these eigenvalues must be negative. The fixed points are 
unstable between the turning points of the state equation, Eq. (19). In Figure [l] we show 



the unstable fixed points as a dashed line in plot (a). Note that from Eq.(20), 



|A±|2^A2 = ^ (26) 
driQ 

and the eigenvalues vanish at the turning points. The linearised analysis thus breaks down 
at the bifurcation points (switching points) 

III. AMPLIFIER GAIN 

Let us now turn to the operation of this device as an amplifier by including the weak 
signal field, 7^ 0. The pump field is chosen to pick out a particular steady state operating 
condition. We then include the signal field in the linearised equations of motion around that 
steady state. In an interaction picture, the equations of motion then become 





(27) 

where 6 = ujg — ujpis the detuning between the signal frequency and the pump. These linear 
equations are easily solved. Here we will be primarily be interested in the signal to noise 
ratio as a function of frequency in the long time limit. To this end we take the Fourier 



transform of Eq. (27) and neglect initial conditions (which decay to zero in the long time 
limit). Define 

poo 

5a{uj) = I e*"*5a(t) (28) 
00 
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The equations of motion then imply that 

5a(uj) \ I ies6(uj — 6) , , , 

\ [ =g{uj) '\ \\ (29) 

where the gain metrix g{uj) is defined by 

g{io) = {M + iujl)-^ (30) 
with / the identity matrix in two dimensions. Then 

r ~o o ^1-1 [ -2-f + Ai(uj + A) AiG \ 

L ^ \^ -AiG* -27 + 4i(tu-A) J 

A = A + Axn-Q. Note that for x = 0, Q is diagonal corresponding to a driven and damped 
linear oscillator. 



The delta function frequency dependance on the right hand side of Eq. (29) is of course 
a consequence of the assumed monochromatic harmonic time dependance of the signal field. 
In general we would have a spectrum for the signal e{uj) in place of e6{uj + 6). In the case 
considered here the signal is harmonic and oscillating at frequency ±5 (in the interaction 
picture). Explicitly, the long time solution in the time domain is 

Sa{t) = tesgii{6)e-"'' - tesg;,{6)e''' (32) 

The appropriate signal quantity is the average of the linearised variable 6X0(1:) = Xg(t) — 
{Xg)o. To be specific we take the local oscillator frequency to be ulo = ujp — 6. The average 
signal then is most conveniently represented in terms of the Fourier amplitude of the mean, 



6x0{uj) = / dte'^^'iSXeit)) (33) 
J —00 

= 6a{u- 6)e''^ + 6a* {u + 6)e-'^ (34) 

where 6a{uj) is the Fourier transform of 6a{t) = {6a{t)) with 6a{t) = ai{t) — ao- The 
homodyne signal is at DC in this picture so we evaluate this at = 0, we find that 

6xg{0) = -te,{[g{6)]ue'' - [g{-6)hie-''} (35) 



where [^(c(j)].jj is a matrix element of g{uj). Using Eq. (31) this becomes 



6xe{0) = —~ ^^^^^^ + c.c. (36) 

^ ^ 4(A2 - ICP) + (7 - 2z5)2 ^ ^ 



where c.c stands for complex conjugate. If we write = noe'^^'^° this becomes 

_ Ae.xnp [jX' - 5') cosjO + 20o) - 7^ sin(g + 20o)] 



where 



A' = Y + (A + 6xTio)(A + 2xno) (38) 



Note that this quadrature amphtude goes to zero when the nonhnearity is turned off. We 
can thus say that this field component is generated by the nonlinear response of the cavity. 
If we now define 

we see that 

Thus for the appropriate choice of local oscillator phase at each bias point, 9 = v — 20o, the 
maximum gain can be defined as 

g[d,no] = , (41) 
v/(A2 - 52)2 + ^252 ^ ^ 

This function is peaked along the lines A = ±5. In Figure |2] we plot the gain versus the 
detuning and no for this optimal choice of phase and with = 7 = 1. Note that for 5 = 
the gain is inversely proportional to square of the eigenvalues of the linearised dynamics. As 
these eigenvalues approach zero at the critical points the gain diverges at this point signaling 
a break down of the linear approximation. A full nonlinear analysis would be required to 
correctly compute the response near the switching points. Determining the conditions for 
maximum gain is not sufficient to specify the operating conditions as the noise might also be 
expected to peak around the bifurcation points. In the next section we turn to a linearised 
noise analysis to determine the signal to noise ratio for the model. 

For completeness we now give the signal when the local oscillator is tuned to ujlo = ujp+S. 
Proceeding as above we find the DC component of the average signal is now given by 

SxeiO) = ^esiQni6)e'' - Q22i-6)e~'') (42) 

This quadrature does not depend on the mean field in the cavity, ao. It thus carries the 
empty cavity response. In fact if we set the pump to zero, turning of the nonhnearity, this 
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(a) Lower 



(a) Upper 




(b) Lower 





(b) Upper 




FIG. 2: Plots of g{5,nQ) versus detuning and no, for the response curves in figure 1 (a) and (b). 
We have separated out the the results for the lower and upper branches in order to avoid the region 
around the critical points where the linearisation breaks down. 



quadrature is given by 



2e. 



52 



(43) 



which is the empty cavity response to the driving field. 



IV. NOISE ANALYSIS 



In order to compute the stationary noise power spectrum for the measured signal we 
again face the problem of not having a specific model of the transducer. However we can 
proceed in general terms by assuming that the transducer is described by a bosonic field 
and the nanomechanical resonator couples the input to this field, amit) to the output aoutif) 
as discussed in the introduction. At high frequencies and for weak coupling the relation 
between the transducer fields and the amplitude of the mechanical resonator is given by 
Eq.g. 

The linearised quantum stochastic differential equations are then 




M 



6a 
5at 



+ 77 



ain{t) 



(44) 
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where we assume that the input transducer field is in a coherent state with vacuum noise 
only, 

{a^n) = {es/Vl)e-''' (45) 
{a,nit),alit')) = 6it-t') (46) 

The complex amplitude of the signal field is coherent so it does not contribute to the noise. In 
what follows we assume that this amplitude has been subtracted off the detected transducer 
signal for the purposes of computing the noise. We thus drop it in the rest of this calculation 
and take the input transducer field to be in the vacuum state. 

Drawing again on the quantum optical analogy with homodyne detection we will simply 
define a stationary noise power spectrum as[ 



/oo 
: 5Xr*(t),5Xr*(0) :) (47) 
-oo 

where the colons denote normal ordering. In effect this means that a shot noise floor has 
been subtracted off the noise power spectrum. In terms of the Fourier component operators 
thus may be written as 



Sr\uj)^ / dJ{: 5Xr\u;),5Xr\u;') :) (48) 

J — oo 

where 

SXr\^) = S^outiu; - 6)e'' + 6dl,{u + 6)e-'' (49) 

Note that when eg = 0, 5 = 0, Eq.([47|) is simply the quadrature squeezing spectrum computed 
in [T7j. However here we are interested in the noise at the signal frequency, i.e. Us- 

We can relate the output quadrature operators of the transducer to the quadrature opera- 



tors of the nanomechanical resonator by taking the Fourier transform of Eq. (44 ) . Combining 
this with the relation in Eq.([9]) we can compute the noise power spectrum of the detected 
signal at the signal carrier frequency 6. This corresponds to evaluating the output spectrum 
at u = 0. 

Expanding the moments in terms of the Fourier amplitudes of the fleld fluctuation oper- 
ators, and applying the appropriate commutation relations, the spectrum can be expressed 
in terms of the gain matrix elements, 

+g2i{uj + 5)g^2{-^ -5) + e-^''g2i{uj + (5)(e?22(-^ - 5) + 1/7)) (50) 
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Evaluating at c<j = gives the noise power in the homodyne signal as a function of signal 
detuning, 



sr(o,5) 



27'|G'|2 + 7 (-iGe^'%{^ - zAf + 6^ + \G\^) + c.c 



(51) 



(A2 - 52)2 + ^2^2 

We can recover the known results for the squeezing in this model [17] when there is no 



injected signal by setting 5 = in Eq. (50) and setting the phase of the local oscillator, 9 



so that the on-resonance noise {uj = 0) is a minimum. 

Here we choose the particular local oscillator phase that was used to fix the gain of the 
average signal, 9 = v — 20O; then 

noutic x\ ^ [8x^^o7 + 4x^o(5^ - + 72/2) sin(z/ - 0o) - 4xno7(A + Axno) cos(z/ - 0o)] 
Sg {0,0} = 



(A2 _ 52)2 + 



(52) 



The phase angles are given by Eq.(24) and Eq.(39). The output noise power spectrum as a 



function of the signal frequency is plotted, in Figure [3j In practice, other phase choices for 
the local oscillator might be preferred; for example one might choose the phase to find the 
quadrature with minimum noise. 

(a) Lower (3) ^PPer 




FIG. 3: Plots of the noise spectrum versus detuning between signal and pump fields, and no, for 
the bias condition shown in curve Figure 1 (a) and (b). We have separated out the the results 
for the lower and upper branches in order to avoid the region around the critical points where the 
linearisation breaks down. 
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We can find an expression for tfie minimum detectable force. First define tlie square of 
tlie square of the signal-to-noise-ratio (SNR), defined as 

Note that in this definition we have added the vacuum noise level to the normally ordered 
variance of the output variable to get the total noise level. It is clear that the gain is peaked 
on 6 = ±A so we set this condition and find the minimum value of for which the SNR is 
unity. 

To evaluate performance we can compare the response of the signal to the nonlinearity 
in the cavity to the response when the nonlinearity is turned off, i.e. the cavity is empty. 
In this case, for our choice of signal phase, the response is given by the quadrature phase 



in Eq. (43) found at dc for a local oscillator tuned to tUp + S. The noise in this quantity is 
simply vacuum noise, so the minimum detectable force, given by setting the signal to noise 
ratio to unity is, 

1 /-Y^ \ 



This is clearly a minimum on resonance as expected, so that the minimum detectable force, 
in units of cavity line width is e™"/7 = 0.25. 

With the nonlinearity turned on, we consider the quadrature phase signal generated by 



the nonlinearity itself, Eq. (37). There are many parameters that one might want to optimise 
for maximum signal to noise ratio. This is made difficult by the fact that the phase of the 
mean amplitude and the phase of the squeezed quadrature are not necessarily the same 
in this model. However for the choice of local oscillator phase we have used here, we see 
that the gain is maximised along the curves 6 = X. Furthermore we expect quantum noise 
reduction inherent in this model to be good close to the instability at A ^ 0. In this limit 
the minimum detectable signal squared is given by 

"'"^^'^ +2xRe(-^^]] (55) 




In Figure |4| we plot e™" versus uq. The dotted line shows the best that can be achieved 
with the nonlinearity turned off. We clearly see a region of operation, close to the switching 
point, in which the minimum detectable force is less than the empty cavity case. 
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FIG. 4: A plot of the minimum detectable force versus no for the mono stabihty case. The dashed 
line represents the minimum detectable force for an empty cavity, a standard quantum limit for 
this model. 

V. CONCLUSION 

In this paper we have proposed a model for the quantum noise of a driven nonlinear 
nanomechanical resonator. We have computed the small signal gain, and signal to noise 
ratio, for a driven and damped nonlinear nanomechanical system at zero temperature. As 
expected the noise is greatest when the system is pumped near the switching point between 
the upper and lower fixed point. Considered as a function of the detuning between the pump 
frequency and the small signal frequency, we see that there are optimal detunings where the 
signal to noise ratio is large. For an appropriate choice of local oscillator phase reference, 
the gain and signal to noise ratio are maximised when the detuning between the pump and 
signal field is equal to the magnitude of the eigenvalues of the linearised dynamics near the 
fixed points. We show that there is a region of operation in which the minimum detectable 
force is less than it would be for a linear cavity. 
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